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We present the three-loop result for the soft anomalous dimension governing long-distance sin¬ 
gularities of multi-leg gauge-theory scattering amplitudes of massless partons. We compute all 
contributing webs involving semi-infinite Wilson lines at three loops and obtain the complete three- 
loop correction to the dipole formula. We find that non-dipole corrections appear already for three 
coloured partons, where the correction is a constant without kinematic dependence. Kinematic de¬ 
pendence appears only through conformally-invariant cross ratios for four coloured partons or more, 
and the result can be expressed in terms of single-valued harmonic polylogarithms of weight five. 

While the non-dipole three-loop term does not vanish in two-particle collinear limits, its contribu¬ 
tion to the splitting amplitude anomalous dimension reduces to a constant, and it only depends on 
the colour charges of the collinear pair, thereby preserving strict collinear factorization properties. 

Finally we verify that our result is consistent with expectations from the Regge limit. 


Infrared (long-distance) singularities are a salient fea¬ 
ture of gauge-theory scattering amplitudes, and a de¬ 
tailed understanding of their structure and how they can¬ 
cel in measurable cross sections is key to precision col¬ 
lider physics. For this reason, there has been a contin¬ 
uous theoretical interest in the factorization and expo¬ 
nentiation properties of the singularities, and their use 
for resummation of large logarithmic corrections, start¬ 
ing from the analysis of the form factor in the early 
days [iHini through to many recent studies of multi-leg 
amplitudes of both massless [3 HUES] and massive par- 
tons P51 - I53 at the multi-loop level, and the formulation 
of the non-Abelian exponentiation theorem in the multi¬ 
leg case dOHll]. 

The focus of this paper will be the infrared (IR) struc¬ 
ture of a scattering amplitude for n massless partons. 
More precisely, if the external legs have momenta Pi, 
i = l..n, with = 0, long distance singularities (both 
soft and collinear) can be factorized as follows 

({Pi} ; ^s) ({Pi} ; ^s) ({Pi} ; T: ^s) ; (1) 

where p is a factorization scale, as = as{p'^) is the 
renormalised U-dimensional running coupling. Pin is a 
finite hard scattering function, and Zn is an operator 
in colour space that collects all IR singularities in the 
form of poles in the dimensional regularization param¬ 
eter e = (4 — Z?)/2. The IR singularities contained in 
Zn have their origin in loop momenta becoming either 
soft or collinear to any of the scattered partons (see e.g. 
Ref. IHj)- Collinear singularities depend on the spin and 
momentum of that particle, and decouple from the rest of 
the process. In contrast, soft (non-collinear) singularities 
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are independent of the spin, but they depend on the rel¬ 
ative directions of motion and the colour degrees of free¬ 
dom of the scattered particles. Hence, soft singularities 
are sensitive to the colour flow in the entire process, and 
their structure is a priori rather complex. Nevertheless, 
they are significantly simpler than finite contributions to 
the amplitude. They can be computed by considering 
correlators of products of Wilson-line operators emanat¬ 
ing from the hard interaction, following the classical tra¬ 
jectory of the scattered particles and carrying the same 
colour charge. 

Specifically, can be obtained as a solution of a 
renormalization-group equation as 



where r„ is the so-called soft anomalous dimension ma¬ 
trix for multi-leg scattering, and P stands for path¬ 
ordering of the matrices according to the order of scales A. 
We stress that r„ itself is finite, and IR singularities are 
generated in Eq. (§ owing to the fact that r„ depends on 
the U-dimensional coupling, which is integrated over the 
scale down to zero momentum. The functional form of r„ 
is highly constrained, and owing to factorization and the 
rescaling symmetry of the Wilson line velocities [T5I120j . 
through three loops it must take the form 

r„ (fe}, A) = rjl'p- ({ft}, A) + A„ ({ftftz}), (3) 

with 

rjl'P- {{P^} , A) = - {as) X log (T, • T, 

Kj ^ ^ 

n 

+ '^ljAas) , ( 4 ) 

i=l 

where —Sij = 2\pi ■ pj \ with Ay = 1 if partons 

i and j both belong to either the initial or the final 
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FIG. 1: Representative 3-loop connected webs contributing to the soft anomalous dimension with 4 coloured lines. 


state and Xij = 0 otherwise; are colour generators 
in the representation of parton i, acting on the colour in¬ 
dices of the amplitude as described in Ref. [n]; 7ic( 
is the universal cusp anomalous dimension [3 |46l |47j, 
with the quadratic Casimir of the appropriate represen¬ 
tation scaled out (Casimir scaling of the cusp anomalous 
dimension holds through three loops |46) ; it may be bro¬ 
ken by quartic Casimirs starting at four loops); jj- are 
the anomalous dimensions of the fields associated with 
external particles, which govern hard collinear singular¬ 
ities, currently known up to three loops [2S1|3S]. Equa¬ 
tion Q is known as the dipole formula, and captures the 
entirety of the soft anomalous dimension matrix up to 
two loops. According to the non-Abelian exponentiation 
theorem [ 33 ] the colour factors in A„ must all correspond 
to connected graphs as shown in Fig. [l] Tripole cor¬ 
rections correlating three partons, with colour factors of 
the form i/“*''^T“T^T^, which could appear starting from 
two loops, are not present in the soft anomalous dimen¬ 
sion at any order because the corresponding kinematic 
dependence on the three momenta is bound to violate 
the rescaling symmetry constraints [T5IBD] . While a con¬ 
stant correction proportional to is excluded 

by Bose symmetry, kinematic-independent corrections in¬ 
volving three lines of the form pbejcde {T“,Tf} 

(last two diagrams in Fig. are admissible and we will 
see that they do indeed appear. The first admissible 
corrections involving kinematic dependence in Eq. 
are then quadrupoles, because four momenta can form 
conformally-invariant cross ratios, 


Pijkl — 


{-Sij){-Skl) 
{-Sik){-Sjl) ’ 


(5) 


which are invariant under a rescaling of any of the mo¬ 
menta. Since diagrams with four colour generators con¬ 
tribute for the first time at three loops, this is the first 
order at which contributions to A„ in Eq. @ may ap¬ 
pear, 


oo ^ 

A„ {{p^jki}) = E (S) {{p^m}) ■ (6) 

Three-loop graphs can connect at most four lines, and 
so the general form of the three-loop correction is com¬ 
pletely determined by the four-parton case and can be 


written as 

({Pijkl}) = 16/af,e/cde| 


E 

l<i<j<k<l<n 


T drjibrriC rr\d '77/ \ 

i jk1 -r {PikjU Piljkj 
■T‘lT\T<^T<}F{p,,kuPuk,) 
T<lT\T^TiT{p,,lk,p^klJ) 


-^E E {T^T^jT^T^}, 


i—1l<j<k<n 

j,k^i 


(7) 


where C is a constant and is a function of two 
conformally-invariant cross ratios. Both C and are 
independent of the colour degrees of freedom. Moreover, 
Eq. Q is the most general three-loop ansatz consistent 
with Bose and rescaling symmetry, so C and T are inde¬ 
pendent of the number of legs n. Note that the terms in 
this sum are not all independent, because of the antisym¬ 
metry of the structure constants and the Jacobi identity. 
A^^^ is independent of the details of the underlying the¬ 
ory and completely determined by soft gluon interactions. 
In particular, this implies that An ^ is the same in QCD 
and in A/" = 4 Super Yang-Mills, and it is therefore ex¬ 
pected to be a pure polylogarithmic function of weight 
five. Its functional form has been constrained by consid¬ 
ering collinear limits and the Regge limit [I81426] . but 
it has so far remained unclear whether three-loop correc¬ 
tions to the dipole formula are present. The purpose of 
the present paper is to compute An . We will present 
its complete functional form, hence determining soft sin¬ 
gularities of any massless multi-leg amplitude at three 
loops. Since C and can be extracted from A 4 , we 
restrict our computation to the case n = 4. Before pre¬ 
senting the final result, we give a brief summary of the 
computation. A complete account of the computation 
will be presented in a forthcoming publication |49j . 

We set up the calculation of the soft anomalous dimen¬ 
sion through the renormalization of a product of semi¬ 
infinite Wilson lines with four-velocities Pk, with /3^ ^ 0. 
By considering non-lighlike lines we avoid collinear sin¬ 
gularities, and obtain kinematic dependence via cusp an¬ 
gles = 2Pi-Pj /We eventually extract A^^^ for 
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FIG. 2: Representative non-connected three-loop diagrams of webs which contribute to the quadrupole term A 4 


(3) 


massless scattering by considering the asymptotic light¬ 
like limit /3^ —^ 0 , where the kinematic dependence re¬ 
duced to conformally-invariant cross ratios as in Eq. ([^. 

In organising the calculation we use of the non-Abelian 
exponentiation theorem and we only compute webs. A 
web can be either an individual connected diagram, as 
in Fig. or a set of non-connected diagrams which are 
related by permuting the order of gluon attachments to 
the Wilson lines [^01444) : representative diagrams from 
such webs are shown in Fig. [21 In either of these cases, the 
contribution to A 4 ^ is associated with fully connected 
colour factors, the classification of which was done in 
Ref. [31] . The sum of all two-line diagrams can be written 
in the form 

G2(1, 2) = dipole - {TJ, T^} {T^, T“} 2), 

where ‘dipole’ stands for a term with a colour factor pro¬ 
portional to Ti •T 2 , which contributes to FJ^'p-. The com¬ 
ponent involving four generators via anti-commutators is 
relevant for the calculation of A„; its kinematic depen¬ 
dence is contained in H 2 { 1 , 2) = i/ 2 ( 2 ,1). Similarly, the 
sum of all three-line diagrams takes the form 


r„ is an operator in colour space that acts on the hard 
amplitude, which is a colour singlet and must therefore 
satisfy [T^ 



H„ = 0 . 


(9) 


Hence, when computing A 4 ^ one may systematically 
eliminate T 4 in favour of Tj, 1 < i < 3, thereby reducing 
all four-line colour factors such as pbefcde^arYbrj.crj.d 

to three-line ones such as |x“, Tf} T^TJ:. As a 

consequence, colour conservation relates sets of diagrams 
that connect a different numbers of Wilson lines. Sum¬ 
ming over all subsets of two and three lines out of four 
and using colour conservation, we find that the function 
J- and the constant C can be expressed in terms of the 
kinematic functions Hn as follows: 


2 r_ 

= H 4 {i,j,k,l) - - H 3 {i,j,k) (10) 

- HsiiJJ) - k) + H 3 {j,iJ) + H 3 [k,{iJ}] 

- H3{k,j,l) - H3{l,i,k) + H3{l,j,k) , 


G3(l,2,3) = r''V^‘'^ ^ {T“,Tf}T5nH3(*,J,fc), 

(ij,fc)G(1.2.3) 

j<k 


with iJ 3 (z,j, fc) = H 3 {i,k,j). The kinematic functions 
H 2 and H 3 are necessarily polynomials in log(— 7 ^) for 
lightlike kinematics where 7^- —>■ —00. Note that in the 
expression for G 3 (l, 2 , 3 ) we omitted the tripole term, 
proportional to which vanishes in this kine¬ 

matic limit. Finally, three-loop webs connecting four 
lines can be cast into the form 


G4(1, 2,3,4) = T“T*T^Tf r‘'^iJ4(l, 2,3,4) -k 


/““/“®Ff 4 (l, 3,2,4) -k ^^ 4 ( 1 ,4, 2, 3) 


( 8 ) 


where kinematic function H 4 satisfies iF4(l,2,3,4) = 
— 114 ( 2 , 1 , 3 , 4 ) = 114 ( 3 , 4 , 1 , 2 ); this function depends on 
logarithms of cups angles as well as on non-trivial func¬ 
tions of cross ratios as in Eq. (j^. 

Another important element in organising the calcula¬ 
tion is colour conservation. The anomalous dimension 




H3{i,j,k) + H3{j,k,i) + H3{k,j,i] 


( 11 ) 


where H 3 {i,j,k) = H 3 {i,j,k) + H 2 ii,j) + H 2 {i,k). The 
above equations put strong constraints on the kinematic 
functions 1 !„: the function J- depends on conformal cross 
ratios, while the functions Hn on the right-hand side of 
Eq. (10 1 depend of logarithms of cusp angles; these must 


therefore conspire to combine into logarithms of confor¬ 
mal cross ratios. In addition, G is a constant, and so the 
kinematic dependence of the functions Hn must cancel in 
the sum in Eq. (11). Our computation satisfies all these 


constraints, which provides a strong check of the result. 

The calculation of the individual graphs is rather 
lengthy, and we will only describe the main steps, defer¬ 
ring a detailed exposition to a dedicated publication [49] . 
We set up the calculation in configuration space, with 
four non-lightlike Wilson lines with four-velocities Pk- 
The position of the three- and four-gluon vertices off the 
Wilson lines are integrated over in D = 4—2e dimensions. 
Following Ref. [351 03], we introduce an infrared reg¬ 
ulator which suppresses exponentially contributions far 











4 


along the Wilson lines. This is necessary to capture the 
ultraviolet singularity associated with the vertex where 
the Wilson lines meet. Upon performing the integral 
over the overall scale, we extract an overall 1/e ultra¬ 
violet pole, and the contribution to the soft anomalous 
dimension is the coefficient of that pole, which is finite 
for each of the diagrams in Fig. (they have no subdi¬ 
vergences) and can be evaluated in U = 4 dimensions. 
Next, we observe that the integrals over the positions of 
the three- and four-gluon vertices give rise to one- and 
two-loop off-shell four-point functions, for which we de¬ 
rive a multifold Mellin-Barnes (MB) representation. Af¬ 
ter integration over the position of the gluon emission 
vertex along the Wilson lines, we obtain a MB represen¬ 
tation of each of the connected graphs for the general 
non-lightlike case, depending on the velocities through 
the cusp angles 7^. In order to proceed, we use standard 
techniques [50] to perform a simultaneous asymptotic ex¬ 
pansion for jij —>■ —00 corresponding to the lightlike 
limit, where we neglect any term suppressed by powers 
of l/7ij- After this procedure, we obtain a collection of 
lower-dimensional MB integrals. The remaining MB in¬ 
tegrals are then converted into parametric integrals using 
the methods of Ref. EH, which can be performed using 
modern integration techniques [52] . The sum over all 
connected graphs is expressible as a linear combination 
of products of logarithms of cusp angles 7ij and single¬ 
valued harmonic polylogarithms [551154] with arguments 
Zijki and Zijki, related to the conformally-invariant cross 
ratios (U by 

Z^ijkl Z^ijkl — Pijkl and (1 Zijki) (1 Z^ijkl) — Pilkj ■ (I^) 


We observe that individual graphs are not pure functions, 
but they involve pure functions of weight five multiplied 
by rational functions in Zijki and Zijki- These rational 
functions cancel in the sum over all connected graphs, 
leaving behind a pure function of weight five, in agree¬ 
ment with the expectation that scattering amplitudes in 
A/” = 4 Super Yang-Mills have a uniform maximal weight. 
Moreover, mixed weight terms do appear in two-line and 
three-line webs, but cancel out in the sum. 

Adding up all the contributions, we find the following 
results for the function and the constant C, 


^)Pijkl^ Pilkj) — -^(f ^ijkl) F{,^ijkl) i 

C* = Cs + 2C2 Ca ) 


(13) 


with 


F{z) — £ioioi(z) + 2C2 [Foqi{z) + Ciqq{z)] , (14) 


where Cwiz) are Brown’s single-valued harmonic poly¬ 
logarithms (SVHPLs) [S3| (see also Ref. [SS]). Note that 
we kept implicit the dependence of these functions on 
z. SVHPLs can be expressed in terms of ordinary har¬ 
monic polylogarithms (HPLs) [54] in z and z. The result 
for F in terms of HPLs is attached in computer-readable 
format to this paper. 


Let us now briefly discuss the main features of the 
final result. First, we note that while F{z) is defined 
everywhere in the physical parameter space, it is only 
single-valued in the part of the Euclidean region (the re¬ 
gion where all invariants are negative) where z and z are 
complex conjugate to each other. Single-valuedness en¬ 
sures that has the correct branch cut structure of a 
physical scattering amplitude [5S1[SS]: it is possible to an¬ 
alytically continue the function to the entire Euclidean 
region while the function remains real throughout [Si- 
Next note that if one considers F{z) as a function of two 
independent variables z and z (not a complex conjugate 
pair) this function has branch points for z and z at 0, 1 
and oo. Crossing symmetry, i.e., crossing some momenta 
from the final to the initial state, is realized in a very 
simple way by taking monodromies around these points. 

Next, let us discuss the symmetries of the final an¬ 
swer for the three-loop corrections to the soft anomalous 
dimension. In the four-line case, Bose symmetry is re¬ 
alised on the cross ratios by the action of the group S 3 
which keeps the momentum pi fixed and permutes the 
remaining three momenta. This group naturally acts 
on the space of SVHPLs by change of arguments gen¬ 
erated by the transformations (z, z) 1—>■ (1 — z, 1 — z) 
and {z,z) >-)■ (1/z, 1/z), with z = Z1234. We note that 
geometrically this symmetry simply acts by exchanging 
the three singularities at z S {0,1,00}. Moreover, it is 
known that the space of all HPLs, and hence also SVH¬ 
PLs, is closed under the action of this S 3 , giving rise to 
functional equations among HPLs, i.e., relations among 
HPLs with different arguments. As a consequence, it is 
possible to express all the terms in Eq. (141 in terms of 
SVHPLs with argument z. 

Besides the action of this group S 3 , there is a second 
symmetry group Z2 acting on the space of SVHPLs. In¬ 
deed, the definition of (z, z) in Eq. (12) is invariant under 
the exchange z ■(->■ z, and hence the function F(z) must be 
invariant under this transformation, i.e., F{z) must be an 
even function: F(z) = F(z). This symmetry is realised 
on the space of SVHPLs by the operation of reversal of 
words, namely, if w is a word made out of O’s and I’s, and 
w the reversed word, then we have C^iz) = Cj^iz) -I- ..., 
where the dots indicate terms proportional to multiple 
zeta values. Even functions then correspond to ‘palin¬ 
dromic’ words (possibly up to multiple zeta values), and 
indeed we see that Eq. (14) is ‘palindromic’. 

Let us now comment on the momentum conserving 
limit of , which is of particular interest because it 
corresponds to two-to-two massless scattering. In this 
limit we have z = z = S12/S13 = —s/{s + t). It follows 
that for two-to-two massless scattering F{z) can be ex¬ 
pressed entirely in terms of HPLs with indices 0 and —1 
depending on s/t, in agreement with known results for 
on-shell three-loop four-point integrals [33 EH] E3 • Fur¬ 
thermore, specialising to the Regge limit and expanding 
Eq. ([14]) at large s/{—t) [53] we find no In^ {s/{—t)) for 
any p > 0: A4 ' simply tends to a constant in this limit. 
This is entirely consistent with the behaviour of a two-to- 
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two scattering amplitude in the Regge limit EIIMIISS]; 
indeed, the dipole formula alone is consistent with pre¬ 
dictions from the Regge limit through next-to-next-to- 
leading logarithms at 0 {al) [50] , 

Finally, let us comment on the behaviour of in the 
limit where two final-state partons become collinear. A 
well-known property of an n-parton scattering amplitude 
is that the limit where any two coloured partons become 
collinear can be related to an (n — l)-parton amplitude: 


MniPl,P2,{Pj}) ^ Sp(pi,p2) A^n-l , (15) 


where P = Pi + P 2 : and Pj are the momenta of the 
(n — 2) non-collinear partons. The splitting amplitude 
Sp {pi,P 2 ) is an operator in colour space which captures 
the singular terms for —)■ 0. All elements in Eq. (15) 
have infrared singularities, and these must clearly be re¬ 
lated. Furthermore, Sp is expected to only depend on the 
quantum numbers of the collinear pair |61) to all orders 
in perturbation theory, and hence also its soft anomalous 
dimension. 


hsp = (r„ - + Asp, (i6) 


must be independent of the momenta and colour degrees 
of freedom of the rest of the process. This property is 
automatically satisfied for the dipole formula, but it is 
highly non-trivial for it to persist when quadrupole cor¬ 
rections are present, as these might introduce correlations 
between the collinear pair and the rest of the process. In 
Refs. [T9l|22] this property was used to constrain A„, but 
this was done under the assumption that C in Eq. (1^ 
vanishes. Given our result for A„ , we may now com¬ 
pute the non-dipole correction to the splitting amplitude 
at three loops: 


Ag= (A(3 )-A®i) 


1||2 


—24 (^5 -I- 2 C 2 C 3 ) 


X 


fabejcde ^rj.a ^ ^ 1^2 


(17) 


('31 

We note that Agp only depends on the colour degrees 
of freedom of the collinear pair, and is entirely indepen¬ 
dent of the kinematics, and hence fully consistent with 
general expectations EH EH- We emphasise that Agp 
is independent of the value of n that was used to com¬ 
pute it, which is remarkable. Indeed, the fact that the 
difference in Eq. 0 is independent of n requires intri¬ 
cate relations between different sets of diagrams and thus 
provides a highly non-trivial check of the calculation. 


To conclude, we have computed all connected graphs 
contributing to the soft anomalous dimension in multi- 
parton scattering and determined the first correction go¬ 
ing beyond the dipole formula. We find that such cor¬ 
rections appear at three-loops already for three coloured 
partons, but they only involve kinematic dependence 
in amplitudes with least four coloured partons, when 
conformally-invariant cross rations can be formed. The 
final result is remarkably simple: it is expressed in 
terms of single-valued harmonic polylogarithms of uni¬ 
form weight five. Finally, we recover the expected be¬ 
haviour of amplitudes in both the Regge limit and in 
two-particle collinear limits, and make further concrete 
predictions in both these limits. 
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